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A  Amplitude  of  transverse  oscillations 

a  Mass  ratio,  (see  C(|.  11-a)) 

Drag  coefficient,  see  E(|.  (1) 
d 

Another  drag  coefficient,  see  Eqs.  (2)  and  {'!) 

C^ll^  Another  drag  coefficient,  see  E(|s.  (6)  and  (7) 

C  Steady-I'low  drag  coefficient 

ds 

Cjj  Mean  in-line  drag  coefficient 
Lift  coefficient 

f.iij  An  inertia  coefficient,  see  E(|.  (1) 

Sill  inertia  coefficient.  See  E(|S.  (2)  and  (r>) 

An  inertia  coefficient,  see  E(|s.  (6)  and  (8) 

D  Cylinder  diameter 

f  A  frequency 

f^  Cylinder  oscillation  frequency 

f^  Natural  freciuency  of  an  elastically-mounted  cylinder 

f*  Vortex  shedding  frequency  for  a  cylinder  at  rest 

L  Length  of  the  test  cylinder 

III  Mass  of  tlie  oscillating  cylinder 

Re  Reynolds  numl)er 

S  Strouhal  nuiiiher 

S  Response  parameter,  material  damping  ratio/mass  ratio 

u 

t  Time 

T  I’eriod  of  transverse  oscillations 

u  Velocity 

U  Velocity 
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Miiximuin  velocity  in  a  cycle 

V  Velocity  of  tlie  aiiibieiit  flow 
Reduced  velocity,  VT/U 

X  Uisplaceinent 

\  =  x/D 

=  x/A 

C  Material  damping  ratio  unless  otlieiA'iise  stated 

V  Kinematic  viscosity  of  fluid 

pj,  Density  of  fluid 

Density  of  the  cylinder,  (apparent  density) 

P^  Fluid  densi  ty/appareiiL  structural  tiensity 

w  Circular  frecpiency 

fre<|uency/natural  freguetjcy 

u  Cyl  inder  osci  1  latloti  frequency/natural  frequency 
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1.  irirucmucriofi 

A.  CrilERAL  REMARKS 

The  linear  or  nonlinear  transient  or  sleatly-slaLe  response  of  a  body 
to  forces  actiny  on  it  has  lony  been  a  source  of  interesLiny  and  practi¬ 
cally  applicable  research.  In  fluid  flow,  the  |)roblem  arises  from  the 
fundamentally  unstable  nature  of  the  shear  layers  and  the  feed-liack 
from  tlie  wake  to  the  shear  layers.  The  resultiny  vortex  motion  yives 
rise,  in  addition  to  a  mean  dray  force,  to  a  siynificant  fluctuatiny 
transverse  force  provided  lliat  part  of  the  hotly  contour  lies  itt  the 
r'eyion  of  the  transverse  pressure  yrailient,  i.e.,  downstream  of  the 
separation  iroints.  Evidently,  the  sliape  of  the  afterbody  and  the  pressure 
yradient  are  related  and  an  afterboily  which  jnevents  the  nonliitear  inter¬ 
action  between  the  two  shear  layers  may  eliminate  tire  transverse  pressure 
yratlieiit  without  prevent  iny  the  syimietric  vortex  foi'iiiation. 

Experiments  have  shown  tliat  elastic  bodies  or  elastically-mounted 
riyiil  bodies  with  a  suitable  afterbody  subjected  to  steady  or  time- 
dejiendent  flow/  can  develop  self-excited  oscillations  with  limit-cycle 
Itehavior  and  that  eneryy  may  be  transferred  from  the  fluid  to  tlie  Iroily 
when  the  vortex-sheddiny  freciuency  brackets  tlie  natural  frerpiency  of 
oscillation  of  the  body.  Thus,  not  only  the  shape  of  the  body  but  also 
its  material  characteristics  (e.y.,  density,  elasticity,  internal 
friction,  tension,  type  of  mountiny,  surface  condition,  yaw  anyle, 
proximity  of  other  bodies  or  free  surface,  etc.)  play  siynificant  roles 
in  the  overall  mechanism.  Some  of  these  variables  controls  the  spanwise 
colierence  of  the  vortices  and  hence  the  inteyrated  transverse  force  alony 
the  body;  some  controls  the  natural  fre<|uency  of  the  body  and  the  mode 

y 


of  its  vibration;  and  su.iie  clianyes  tl.e  internal  damping  of  the  body  and 
the  work  done  by  its  sun  oinulimjs.  Evidently,  means  may  be  discovered 
through  painstakingly-di f ficul t  experiments  to  leduce  or  destroy  the 
spanwise  coherence  so  as  to  minimize  the  exciting  force,  to  increase  the 
internal  damping  so  as  to  minimize  the  net  energy  transfer,  or  to  increase 
the  natural  frequency  (stiffening  of  the  body  or  of  its  supports)  so 
tl>aL  the  vortex-shedding  f rcH|tieucies  are  well  l)elov/  or  aliove  the  natural 
frequency  of  tlie  body.  In  doing  so,  however,  Llie  objective  is  not  only 
the  reduction  or  elimination  of  the  hydroelasLic  oscillations  but  also 
the  minimization  of  an  in-line  drag-force  parameter  such  as  (l^A/l))CJ. 

In  other  words,  tlie  penalty  for  the  reduced  transverse  oscillations 
should  not  be  a  large  in-line  drag  force. 

The  phenomenon  is  relatively  more  complex  since  one  is  concerned 
not  only  with  the  oscillations  of  the  body  and  Die  in-line  and  transverse 
forces  acting  on  it  but  also  with  fatigue  and  the  sound  generated  by  the 
oscillations  and  vortices.  I'urthermore,  in  affecting  the  natural  frequency 
of  the  body  by  additional  tension  or  mass,  the  (|nestion  is  not  merely  the 
deteniiination  of  the  separation  needed  between  the  vortex-shedding  and 
natural  frequencies  but  also  of  the  range  over  whicli  the  excititjg  mechanism 

I 

can  lock  onto  the  frequency  of  tlie  body. 

The  material  damping  which  is  often  combined  or  confused  with  fluid 
damping,  cannot  be  increased  i.j  desired  levels  without  severe  penalties. 

As  to  the  efforts  to  change  the  body  shape,  it  is  obvious  that  a  circular 
cross-section  is  most  desirable  partly  because  of  ease  of  Us 
construction  and  partly  because  of  the  fact  that  often  the  flow  is 
omnidirectional  as  well  as  t ima-ile|)endenL.  It  is  for  the  reasons  cited 
above  that  the  reduction  of  the  integrated  transverse  force  by  a  reduction 
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of  tlie  spanv/ise  coherence  (increasiny  the  force-cancell  iny  effects  of 
phase  sliifts  which  may  be  bronyht  about  by  three-dimensional  effects) 
by  means  of  iielical  strakes,  eccentric  rinys,  uneven  rouyhness,  and 
various  types  of  fai rinys  has  been  the  primary  method  of  vibration 
suppression.  It  should  be  noted  that  such  devices  are  often  referred  to 
as  the  "vortex  suppression"  devices  even  thouyh  they  do  not  actually 
suppress  the  vortex  sheddiny.  Tliey  merely  reduce  the  integrated  lift 
force  through  fihase  shifts. 

It  is  evident  from  tlie  foregoing  remarks  that  there  are  substantial 
yaps  in  the  understanding  of  the  fundamental  mechanisi.is  at  v/ork  in  the 
forcing  of  a  body  l)y  the  forces  acting  on  its  afterbody  and  in  the 
predictioti  of  the  lock-in  range  in  which  resonant  and  pseudo-resonant 
oscillations  occur.  It  is  also  evident  that  the  control  or  prevention 
of  the  oscillations  requires  first  of  all  the  prediction  of  whether  such 
oscillations  will  at  all  occur  and  if  so  for  what  range  of  tiie  controlling 
j)arameters. 

In  view  of  the  foregoing  considerations,  the  present  research  program 
v;as  undertaken  with  two  main  olijectives:  (a)  to  determine  tlje  in-phase  and 
out-of-pliase  components  of  tlie  time-dependent  force  acting  on  a  rigid 
circular  cyl inder  undergoing  forced  transverse  oscillations  in  a  uniform 
stream  and  to  predict  through  the  use  of  these  force  components  the 
dynamic  response  of  an  elastically  mounted  cylinder  for  which  tlie  material 
damping  coefficient  and  tlie  'mass  parameter'  are  known;  and  (b)  to  predict 
through  the  use  of  the  discrete  vortex  model  the  kinematic  and  dynamic 
characteristics  of  the  fluid  motion  about  a  circular  cylinder  undergoing 
harmonic  oscillations.  The  second  objective  viill,  hoiiefully,  help  to 
synthesize  the  information  gathered  during  the  jiast  decade  from  specific 
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practical  vibration  problems  and  from  iiigiily  idealized  experimental  and 
theoretical  investigations. 

This  report  does  not  deal  with  vibration  supi)ression  methods,  the 
effect  of  time-dependence  and  shear  of  the  ambient  flow,  galloping  motion, 
the  effect  of  yaw,  coupled  in-line  and  transverse  oscillations  of  flexible 
cylinders  or  cables,  the  proximity  effect  of  other  tubes  or  bodies,  etc. 
lurthermore,  no  attem|)t  iias  been  made  to  offer  a  chronological  survey  or 
a  state-of-the  art  appraisal  of  the  flow  inducetl  osci  1  latiofis.  Only  those 
v/orks  which  have  a  direct  liearing  on  the  evaluation  and/or  discussion  of 
the  jjresent  results  are  reviewed  in  some  detail  v/herever  appropriate. 

11.  SEl  F-EXClTrO  OSCli.l  AI  lOiiS 

The  net  result  of  alternate  vortex  sliedding  is  an  oscillating  side 

thrust,  upon  a  cylifider  of  suitable  form,  in  a  direction  away  from  the 

last  vortex.  This  side  thrust  or  lift  force  exists  practically  at  all 

Ileynolds  numbers  regardless  of  whether  the  body  is  allowed  to  respond 

dynamically  or  not.  For  a  body  held  at  rest,  there  is  a  finely-tuned, 

relatively  delicate,  interaction  and  lialaiicc  between  various  features  of 

the  flov/.  tJ|)Stream  bouinlary  layers  generate  the  vorticity  which  feeds 

1 

the  vortices.  Downstream  boundary  layers,  interaction  between  the 
oppositely-signed  shed  vortices,  and  turbulence  control  the  dissipation 
mechanism,  the  amount  of  vorticity  each  vortex  receives  or  retains,  the 
mobility  of  the  separation  points,  the  average  vortex  sliedding  frequency, 
etc.  The  fact  that  the  various  components  of  such  a  flow  is  delicately 
interrelated  and  susceptible  to  interference  has  been  amply  demonstrated. 
The  appreciation  of  the  fact  that  the  alternate  vortex  shedding  is  a 
consequence  of  flow  instability,  or  of  the  break«lown  of  the  cause  and 
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effect  syiiuiietry,  is  particularly  important  in  understanding  as  to  why 
the  vortex-induced  oscillation  is  not  a  simj)le  cause-and-effect  phenomenon. 
From  a  fundamental  point  of  view  the  question  is  the  understanding  of 

both  the  modification  of  the  steady  flow  mechanisms  and  the  interdependence 

1 

between  the  resulting  time-dependent  force  that  causes  the  body  to  oscillate 
and  the  oscillation.  From  a  practical  point  of  view  the  question  is  how 
to  prevent  the  body  from  taking  control  of  vortex  shed<ling. 

Numerous  ex|)eriments  have  shov/n  that  (e.g.,  Blevins  1977)  when 
the  vortex  shedding  frec|uency  brackets  the  natural  frequency  of  a  body  v/ith 
a  suitable  afterbody,  the  body  takes  control  of  the  shedding  in  apparent 
violation  of  the  Strouhal  relationship.  Furthermore,  the  spanwise 
correlation  of  the  wake,  the  vortex  slrength,  ami  the  drag  force  increase. 

The  base  pressure,  wake  width,  vortex -shedding  frecpiency  may  increase  or 
decrease  as  one  approaches  the  synclironization  region.  These  are  merely 
some  of  the  s imply-obsei  val)le  or  measurable  flow  features  vdiicl)  in  reality 
reflect  the  changes  in  the  fine-grain  flow  structure  due  to  the  complicated 
interaction  between  the  oscillator  and  the  oscillated.  Thus,  it  may  he 
said  that  the  lock-in  iJhenomenon  over  a  range  of  velocities  owes  its 
occurrence  to  the  precarious  nature  of  the  base  flow.  The  more,  |»recarious 
the  basic  flow  features  are  al)out  the  body,  the  more  susceptible  is  the 
body  to  synchronous  oscillations  in  a  v/ider  velocity  range.  In  other 
words,  a  body's  ability  to  alter  the  flow  features  through  its  oscillations 
or  simply  to  amplify  the  lift  force  is  a  measure  of  its  ability  to 
hydrodynamic  response.  It  is  not  reasonable  to  think  that  the  flow 
characteristics  about  an  oscillating  body  woulil  be  the  same  as  those  about 
a  body  at  rest.  Should  tiiis  have  l)een  the  case,  there  would  have  been 
resonance  at  one  particular  velocity  of  flow  only,  namely,  when  the  Strouhal 
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frequency  is  equal  to  tlie  natural  frequency  of  tlie  body. 

The  interplay  betv/een  the  different  piiysical  quantities  involved  in 
the  phenomenon  can  be  adequately  established  and  a  deeper  insiylit  could 
be  gained  into  the  physics  of  the  floi»  only  through  an  analytical  model. 

To  a  degree,  the  difficulties  in  devising  a  model  stem  from  an  inadequate 
knowledge  of  tl|e  precise  mechanism  of  flow  sei>aration  and  its  consequences 
in  steady  flow  past  a  stationary  body.  To  be  sure,  analytical  models  in 
vdiich  the  e(|uations  of  motion  are  conditioned  in  a  manner  analogous  with 
experimental  tecliniques,  such  as  conditioned  sampling,  for  instance, 
could  also  assist  in  exposing  the  details  of  the  piiysical  piienomena.  On 
the  other  hand,  models  contrived  to  suit  tiie  available  experimental  data 
give  only  an  illusion  of  understanding  of  the  |)hysics  of  the  flow  and 
become  at  best  an  interiiolation  lietv/een  various  experimental  results. 

So  far,  it  has  not  been  possible  to  calculate,  by  appropriate  viscous 
or  inviscid  modeling,  the  development  of  coiled-up  vortex  sheets  springing 
from  two-  or  three-dimensional  separation  lines  on  a  body  of  general  sliape. 
Thus,  isolating  a  small  number  of  flow  features  either  to  describe  tlie  wake 
development  about  a  stationary  boily  or  the  differences  between  the  charac¬ 
teristics  of  stationary  and  oscillating  bo<ly  so  as  to  arrive  at, a 
theoretical  unification  of  the  existing  observations  and  ineasurbments  may 
be  an  illusive  search.  However,  the  end  result  may  be  less  important  than 
the  practical  understanding  that  the  effort  lias  brought.  In  the  main, 
it  seems  that  the  computation  of  tlie  characteristics  of  the  near  v;ake 
region  and  tlie  interaction  of  that  region  witli  an  oscillating  body  will 
continue  to  remain  essentially  empirical  and  the  success  of  numerical 
schemes  v/ill  be  sensitive  to  the  choice  of  i/idependent  parameters.  It 
appears  that  heavy  reliance  on  experimental  investigations  of  vortex- 
induced  oscillations  and  the  evolution  of  essentially  empirical  models 


14 


/ 

to  brid'je  between  sets  of  experimental  results  will  doiainate  tl«e  problctn 
of  fluid-structure  interaction  for  a  lony  time  to  come. 

C.  Discussion  OF  RCLCVAItr  PHYSICAL  PAUAHLTFUS 

Flov/  about  an  oscillating  body  necessarily  concerns  such  parameters 
as  the  Strouhal  number,  damping  coefficient,  natural  frequency,  added  mass, 
etc.  These  parameters  are  variously  interpreted,  evaluated,  or  assumed 
by  various  workers.  Thus,  it  is  necessary  that  they  be  discussed  here 
briefly  so  tliat  tlie  reader  can  assess  the  validity  of  some  of  the 
assumptions  made  in  their  use  in  various  analytical  models. 

a.  Strouiial  Number 

The  separation  process  is  a  time-dependent  phenomenon  with  alternate 
shedding  vortices  creating  a  transverse  pressure  gradient.  Unsteady 
hydrodynamic  loads  arising  from  tliis  pressure  gradient  acting  on  tlie 

I 

afterbody  can  excite  dynamic  structural  response  as  noted  earlier,  body 

natural  frecpiencies  near  the  exciting  frequency  raise  the  spector  of  load 

and  response  enhancement.  Thus,  the  spectral  content  of  the  forcing 

functions  are  important  to  dytiamlc  structural  rospotise  analysis.  If  the 

shedding  process  is  periodic  or  quasi -periodic ,  a  simple  characteristic 

Strouhal  number  S'*  =  f“D/U  is  defined  where  f"  is  the  shedding  frequency 

V  v 

when  the  cylinder  is  at  rest.  The  shedding  process  may  be  random 
(broadband)  over  a  portion  of  the  Reynolds  number  range  for  winch  a 
statistical  response  analysis  is  required  based  on  the  spectral  content. 

Also  required  is  a  measure  of  the  correlation  length  along  the  cylinder. 

If  the  correlation  lengths  are  long  compared  to  the  cylinder  diameter, 
sectional  analyses  are  valid,  and  the  loads  act  in  concert  along  the 
cylinder.  If  short  correlation  lengths  are  present,  the  structural 
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analysis  must  be  three-diineiisional ,  and  tlie  total  load  on  the  body  Is 
reduced.  Evidently,  tlie  characterization  of  the  vortex  shedding  process 
by  a  simple  frequency  is  a  practical  simplification.  Pov/er  spectral 
density  analyses  of  unsteady  cylinder  loads  reveals  that  in  certain 
Reynolds  number  regimes  the  shedding  process  is  periodic  and  can  be 
characterized  by  a  single  freepjency.  At  subcritical  Reynolds  numbers, 
the  energy  containing  frecjnencies  are  confined  to  a  narrov<  band,  and  tlie 
Strouhal  number  is  about  O.'^?  lor  smooth  cylinders.  It  must,  however,  be 
emphasized  that  only  an  average  Stroulial  numlier  may  be  defined  for  Reynolds 
numbers  larger  than  about  20,000.  In  the  critical  Reynolds  number  regime, 
a  broad  band  power  spectral  dcMisity  is  usually  oliserved.  At  higher 
Reynolds  numbers,  the  Stroulial  fretiuency  rises  to  about  0.3  and  the 
shedding  process  is  quasi -periodic.  The  spectral  content  of  the  exciting 
forces  is  particularly  important  for  bodies  which  may  undergo  in-line  and/or 
transverse  oscillations  since  the  vortex  shedding  frequency  locks  onto  the 
frequency  of  the  transverse  oscillations  of  the  cylinder. 

Various  attempts  have  been  made  to  devise  a  universal  Strouhal  number 
which  would  remain  constant  for  differently  shaped  two-  or  three-dimensional 
body  shapes  (e.g.*  Simmons  197/).  Unfortunately ,  these  definitions  suffer 

i 

from  the  obvious  drawback  that  their  definition  requires  tiie  solution  of 
the  wake  formation  problem  first,  or  the  measurement  of  one  or  more  flow 
characteristics.  Thusi  their  importance  lies  not  so  much  in  their  ability 
to  predict  but  rather  to  uncover  the  intricate  relationship  between,  say, 
the  flow  velocity,  vortex  shedding  frequency,  wake  width,  base  pressure, 
the  streaimvise  distance  from  separation  edges  to  the  point  of  minimum 
base  pressure,  or  the  distance  apart  of  the  peaks  of  the  r.m.s.  value  of 
the  velocity  across  the  wake,  etc.  It  appears  that  isolating  a  small 
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lumiber  of  simple  scales  to  describe  the  wake  development  may  not  be 
poss ible. 

As  noted  earlier,  the  Sttouhal  number  appears  to  be  a  function  of 
the  Reynolds  number.  The  fact  that  inviscid  flow  models  such  as  the 
discrete  vortex  analysis  predict  the  Strouhal  number  almost  exactly 
in  conformi ty -wi th  the  cyueriments  in  the  sober itical  Reynolds  number 
ranye  is  rather  surprising. 

The  constancy  of  the  Strouhal  number  over  a  broad  range  of  Reynolds 
numbers  does  not  imply  that  the  base  pressure  remains  constant  and  that 
a  single  vortex  emanates  from  a  separation  line  each  time  a  vortex  is  shed. 

In  reality,  there  is  not  only  a  phase  shift  between  various  sections  of 
the  vortex,  separated  by  a  correlation  length,  but  also  variations  in  both 
tlie  intensity  and  the  frecpiency  of  the  vortex  segments.  The  variation 
of  the  liase  pressure  with  Reynolds  number  in  the  range  where  the  Strouhal 
number  remains  constant  may  be  related  to  the  variation  of  the  correlation 
length  with  tlie  Reynolds  number.  The  net  effect  of  the  spanwise  variations 
of  the  vortex  tube  is  that  the  lift  coefficient  olitained  from  a  pressure 
integration  is  not  necessarily  identical  v/ith  that  obtained  from  the  direct 
total  force  measurements.  Partial  spanwise  correlation  leads  to  variations 

I 

in  both  the  frequency  and  the  amplitude  of  the  lift  force,  the  variation 
of  the  latter  being  more  pronounced  than  that  of  the  former.  The  reasons 
for  these  variations  and  the  lack  of  total  spanwise  correlation  are  not 
quite  clear.  The  end  effects,  wal 1 -boundary  layers,  freestream  turbulence, 
non-uniformity  of  the  velocity  distribution  are  mentioned  quite  often  as 
possible  reasons.  It  is  possible  that  the  flow  along  tlie  cylinder  is 
Just  as  unstable  as  it  is  cross  to  it.  The  instability  of  the  shear  layers 
to  a  band  of  frequencies  and  the  amplification  of  the  disturbances  in 
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response  lo  various  frequencies  alony  the  cylinder  could  lead  to  a  three 
dimensional  flow.  This,  in  turn,  may  yive  rise  Lo  a  complicated  interaction 
between  the  near  wake  and  the  mobile  separation  point  on  a  circular  cylinder. 

It  would,  hov/ever,  be  wrong  to  assume  that  the  mobility  of  the  separation 
points  is  primarily  responsible  for  the  imperfect  spanwise  coherence.  Even 
bodies  such  as.  90  degree  wedges,  s(|uare  cylinders,  etc.,  with  fixed 
separation  lines  (assuming  no  reattachment)  do  not  exhibit  perfect  correlation. 
To  be  sure,  the  variation  of  the  l)ase-pressure  coefficient  for  bodies  with 
mobile  separation  lines  is  greater  than  that  for  bodies  v/ith  fixed  separation 
points  (Roshko  1970). 

Suffice  it  to  note  that  there  is  a  lack  of  coherence  in  the  data  for 
the  lift  coefficient  for  a  cylituler  at  rest  l)ecause  of  the  various  reasons 
discussed  above.  The  recognition  of  the  flow  instabilities  and  the  three- 
dimensional  nature  of  tlie  base  flow  yive  a  greater  appreciation  of  the 
difficulties  encountered  in  assessing  tlie  influence  of  the  dynamic  response 
of  the  body  on  the  entire  flow  structure  and  iri  ilcvising  suitable  models. 

b,  Uampifig  Coefficient 

It  is  a  well-known  fact  that  all  systems  dissipate  energy  through 
various  mechanisms.  The  ability  of  the  system  to  do  so  is  called  damping. 

Its  rate  depends  on  the  internal  friction  of  the  material,  the  support 
conditions,  and  tlie  surrounding  medium.  The  effect  of  the  internal  or 
external  fluid  medium  may  be  to  increase  the  apparent  mass,  fluid  friction, 
flov;  separation,  fluid  sloshing,  or  a  combination  thereof.  From  a 
practical  point  of  view,  it  is  simpler  and  often  more  desirable  to 
determine  the  overall  damping  of  a  structure  in  the  fluid  medium  in  vdiich 
it  may  be  subjected  to  hydroelasLic  oscillations.  This,  however,  may  cause 
several  mathematical  and  conceptual  difficulties.  Firstly,  the  part  of 
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the  (lamping  vdiich  Is  attributed  to  tlie  surrouiuiiiKj  fluid  is  also  part 
of  the  driving  or  resisting  force,  depending  on  the  relative  directions 
of  the  fluid  force  and  the  body  velocity.  Thus  the  inclusion  in  the 
dam[)ing  term  of  a  constant  damiring  coefficient  would  be  wrong  since  the 
fluid-induced  damping  is  not  and  cannot  be  constant.  On  the  other  hand, 
the  partition  of  the  effect  of  the  fluid  motion  partly  as  fluid  damping 
and  partly  as  drivingor  resisting  force  would  make  the  matters  worse 
or  ho|)elessly  complex.  The  inclusion  of  the  entire  effect  of  the  fluid 
motion  in  the  damping  term  v/ould  recpiire  the  assumption  of  a  hicjhiy 
nonlinear  damping  coefficient.  This,  in  turn,  would  give  rise  to 
additional  assumptions  and  perhaps  to  unexplainable  conclusions.  In 
spite  of  these  arguments,  it  has  often  been  assumed  that  the  damping 
coefficient  (and  the  natural  fre(|uency)  may  be  determined  by  plucking 
excitation  of  the  cylinder  in  the  fluid  medium.  Oscillations  created 
in  this  manner  are  of  small  amplitude  and  often  free  from  vortex  shedding, 
if  not  from  flow  separation.  For  A/1)  smaller  than  about  unity,  the 
effect  of  the  fluid  motion  is  to  increase  the  linear  friction  force  and 
the  mass  of  the  body  by  the  added  mass  (nut  necessarily  the  displaced  mass 
(Sarpkaya  1976).  On  this  basis,  it  is  often  assumed  that  frictional 
effects  of  the  fluid  motion  may  be  Included  in  the  overall  daniping 
coefficient  and  the  added  mass  effects  may  be  accoutited  for  by  properly 
decreasing  the  natural  fre(|uency.  When  a  body,  for  example  a  cylinder, 
is  subjected  to  transverse  oscillations  in  a  steady  flow,  the  fluid  forces 
resulting  from  vortex  shedding  (essentially  form  drag)  are  taken  as  the 
effective  driving  fluid  force.  The  apparent  success  of  this  procedure 
owes  its  success  to  the  relatively  small  amplitude  oscillations  observed 
and  partly  to  the  difficulty  of  determining  the  total  force. 


The  use  of  a  similai'  procedure  for  ll»e  added  mass  causes  even  more 
difficult  problems  as  v/ill  be  discussed  later.  There  is  no  reason  in 
the  lavjs  of  fluid  dynamics  that  tlie  added  mass,  or  the  drifted  mass, 
should  be  the  same  for  a  body  osci  1  lat iiiy,  at  its  natural  frequency, 
in  a  fluid  otherwise  at  rest  v/ith  tliat  oscillating  in  a  fluid  in  motion. 

It  is  of  course  recognized  that  it  Is  rather  difficult  to  determine  the 
material  damping  of  a  body  in  vacuum  and  the  added  mass  as  a  function  of 
time  in  a  time-dependent  flov/. 

c.  natural  Fre(|uency  and  Added  Mass 

Natural  frequency  is  perhaps  the  most  im|»oi  tant  parameter  related 
to  the  dynamic  response  of  a  system  since  the  synclironization  takes  place 
at  or  near  to  it.  The  determination  of  the  natural  frequency  of  a  body 
it>  air  (say  by  plucking  excitation)  is  (piite  correct  since  the  added  mass 
due  to  the  displaced  mass  of  air  is  much  smaller  than  the  mass  of  the 
body.  In  water,  however,  the  added  mass  may  be  (piite  large  and  l.he 
apparent  natural  frequency  may  be  significantly  lower.  The  justification 
for  the  use  of  a  natural  frecpiency  measured  by  plucking  excitation  in  a 
fluid  otherv/lse  at  rest  is  necessarily  based  on  the  assumption  that  the 
added  mass  is  approximately  ecpial  to  its  inviscid  flow  value  (not  the  mass 
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of  the  fluid  displaced  by  the  body!)  regardless  of  mode  shape,  vibration 
amplitude,  or  presence  of  vortex  shedding.  Strictly  s|)eaking,  what  is 
added  to  the  mass  of  the  body  is  not  the  "mass  of  fluid  displaced  by  the 
body"  but  rather  the  added  mass  olitained  by  classical  potential  flow 
methods  for  the  particular  direction  of  the  motion  of  the  body.  For 
example,  the  displaced  mass  of  a  thin  rectangular  plate  is  practically 
zero.  Its  added  mass  may  be  equal  to  the  added  mass  of  the  cylinder 
accepting  the  plate  as  its  diameter  or  to  zero  depending  on  whether  the 
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plate  is  oscillating  tionnal  or  parallel  to  its  surface.  It  is  also 
important  to  note  that  the  adtleci-mass  and  the  added-iiiass  moinenl  of  inertia 
coefficient  depend,  in  addiliun  to  the  direction  and  type  of  motion,  on 
tl«e  proximity  of  other  bodies  and  the  free  surface. 

The  acceptance  of  the  assumption  that  the  added  mass  does  not 
depend  on  the  mode  of  vibration,  vibration  amplitude,  or  presence  of 
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vortex  shedding  may  lead  to  unexplainable  predictions.  Ttie  added  mass 
coefficient  obtained  by  plucking  excitation  of  a  cylinder  in  a  fluid 
olherv/ise  at  rest  does  not  prove  that  it  remains  constant  for  all  types 
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of  fluid  motion.  It  rather  proves  that  the  added  mass  calculated  from 
potential  flow  theory  is  essentially  correct  for  the  conditions  approximating 
the  potential  flow. 

Added  mass  is  one  of  the  least  understood  and  most  confused  charac¬ 
teristic  of  time-dependent  flows.  Several  facts  need  to  be  made  clear. 
Firstly,  added  mass  exists  regardless  of  whether  the  fluid  and/or  the 
body  are  accelerating  or  not.  It  is  a  drift  of  mass  resulting  from  the 
"elastica  shaped"  path  of  the  fluid  particles  during  the  motion  of  the 
fluid  about  the  body  (Darv/ifi  1953).  This  mass,  like  all  masses,  reveals 
its  existence  only  when  it  is  subjected  to  an  acceleration.  Thus,  the 
inertial  force  required  to  accelerate  the  fluid  about  a  body  at  rest  is 
a  sum  of  the  force  required  to  accelerate  the  added  mass  and  the  force 
acting  on  the  displaced  mass  of  the  body  due  to  the  pressure  gradient 
causing  the  fluid  to  accelerate.  The  inertial  force  required  to  accelerate 
a  body  in.  a  fluid  otherwise  at  rest  is  a  sum  of  the  force  required  to 
accelerate  the  mass  of  tlie  body  and  the  force  to  accelerate  the  added  mass. 

Secondly,  added  mass  depends  on  tlie  type  of  motion  of  the  body  or 
of  the  fluid  about  tlie  body;  proximity  of  other  bodies,  free-surface,  etc.. 


and  time.  It  is  not  always  possible,  nor  advisable,  to  determine  the 
instantaneous  value  of  the  added  mass.  For  example,  for  a  periodic 
flow  cl»aracterized  by  u  =  d^^sin  2i(t/T,  tl«e  added  mass  coefficient 
cannot  be  properly  defined  at  certain  times. 

Thirdly,  the  added  mass  coefficient  (or  its  time  averaged  value 
over  a  suitable  time  interval)  depeinls  on  the  vortex  shedding.  For  a 
cylinder  undergoing  liarmonic  oscillations  with  a  relative  amplitude  of 
A/I)  nearly  etpial  to  unity  in  a  fluid  otherwise  at  rest,  there  is  no 
vortex  s(iedding  and  the  added  mass  coefficient  is  etpial  to  its  classical 
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value  of  unity.  For  larger  values  of  A/I)  where  there  is  vortex  shedding, 
the  added  mass  coefficient  may  even  be  negative  (sarpkaya  1976). 
Conseipiently ,  it  is  wrong  to  equate  tlie  two  flow  situations  where  in  one 
case  the  cylinder  oscillates  with  amplitudes  of  A/U  smaller  than  unity 
in  a  fluid  otherwise  at  rest  and  the  other  case  where  the  cylinder 
oscillates  at  similar  amplitudes  in  a  transverse  direction  in  steatly 
flow.  Tlie  former  does  not  involve  vortex  slicdding  whereas  tlie  latter 
is  accompanied  by  complex  separation  and  vortex-shedding  phenomena 
even  though  the  amplitude  of  oscillations  is  no  greater.  The  empirical 
models  arrived  at  by  ignoring  these  well  estalilished  facts  require 
variable  empirical  constants  to  simulate  the  particular  set  of  data  and, 
as  noted  before,  become  at  best  an  interpolation  between  various  sets  of 
data.  ‘ 

I 

^  t 

Suffice  it  to  note  that  the  modeling  of  the  separated  flow  about  an 
oscillating  body  requires  first  of  all  the  elucidation  of  the  correct 
physics  and  the  separation  of  the  fluid  forces  acting  on  the  body  from 
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such  quantities  as  the  actual  mass,  material  damping,  and  the  spring 
constant  of  the  body.  Tlie  empirical  models  (i.e.,  van  der  Pol  oscillator 
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n«(Iel  or  its  improved  versions,  v/ake-osci  1  lator  models  based  on  a  bidden 
flow  variable  and  the  integral -momentum  etpiation,  etc.)  are  highly 
conditioned  by  the  available  experimental  data.  Thus,  their  power  of 
prediction  is  not  otily  uncertain  but  also  limited  to  the  range  of  the 
experimental  data  used  in  tiieir  evolution,  it  is  also  unfortunate  tliat 
the  empirical  models  obscure  the  physics  of  the  flow  and  make  it  imi)ossible 
to  establish  a  meaningful  relationsliip  between  certain  observed  or  measured 
features  pf  tl>e  near  wake  (e.g.,  narrowing  or  widening  of  the  wake,  the 
increase  'or  decrease  of  the  wake  establishment  region,  tlie  increase  or 
decrease  of  the  strength  of  the  sited  vortices,  etc.).  As  stated  earlier, 
liowever,  tlie  oscillation  of  a  body  in  a  fluid  in  motion  is  a  furtlier 
complication  added  onto  an  already  complex  and  relatively  unstable  flow 
situation.  Thus,  it  does  not  appear  that  lliere  will  ever  be  an  exact 
solution  or  a  comprehensive  model  capable  of  predicting  the  observed 
characteristics  of  tite  plienomenon.  After  all,  tiiose  wiio  liave  concentrated 
on  semi-empirical  ntodels  may  have  been  tlie  wiser  in  recognizing  the  need 
to  get  to  tiie  designer  before  the  understanding  of  the  problem  arrives. 
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II.  EXI’ClUMENTAL  EqUlPMCNF  AND  PROCEDURES 


The  experiments  v/ere  performed  in  two  separate  water  tunnels. 

The  first  one  was  a  recirculating  water  tunnel  with  a  capacity  of 
approximately  500  gallons.  The  galvanized  test  section,  four  indies 
wide,  eight  indies  iiigli  and  16  indies  long,  was  closed  on  to|)  v/itli  a 
removable  plexiglass  plate  to  eliminate  tlie  free  surface  effects.  A 
small  space  adjacent  to  tlie  side  walls  vias  provided  to  allow  tlie  passage 
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of  tliin  leaf  arms  coniiecLiiiy  tlie  cylinder,  in  the  test  section,  to  the 
driving  liardware  above  I  lie  test  section.  A  low  RI’M,  liigli  capacity, 
1'1-incli-diaineter  disdiarge  centrifugal  pump  was  used  to  circulate  tlie 
fluid  througli  tlie  test  section.  The  velocity  of  the  fluid  was  regulated 
by  a  butterfly  valve  arraiigeiiient  situatcil  downstream  of  the  test  section. 
Velocities  of  0.72  to  1.55  feet-per-second  were  obtained  by  adjusting 
tlie  vanes  of  the  butterfly  valve. 

The  periodic  motion  was  obtained  by  use  of  a  small  ,varial)le  speed, 
electric  motor  and  flywheel  and  pivot  asseiiilily.  Tlie  amplitude  A  of  tlie 
vertical  periodic  motion  was  set  by  adjusting  tlie  radial  position  of  the 
bearing  attached  to  the  flywlieel.  For  tliese  experiments  the  amplitudes 
ranged  from  0.25  to  0.R1  inches.  Tlie  frequency  of  oscillation  was 
regulated  by  adjustment  of  tlie  variable  speed  motor.  The  system  was  . 
designed  and  constructed  to  produce  essentially  sinusoidal  oscillation, 
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free  of  secondary  oscillations. 

The  test  specimens  used  were  0.7  indies  and  1.0  inch  in  diameter 
and  constructed  from  aluminum  tubing.  The  cylinder  vias  held  in  the  test 
section  by  a  yoke  assembly  connected  to  the  pivot  arm  by  a  vertically 
constrained  rod.  The  action  of  tlie  pivot  arm  caused  the  cylinder  to 
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oscillate  vertically  In  the  test  section,  transverse  to  the  flow.  The 
aluminum  yoke  assembly  was  instrumented  with  strain  gages  and  accelerometers 
to  monitor  the  forces  and  the  acceleration  felt  by  the  cylinder. 

The  mean  velocity  of  the  fluid  was  determined  by  use  of  a  pitot  tube 
and  a  hot-film  probe  installed  in  the  test  section  upstream  of  the  test 
cylinder.  The  velocity  was  continuously  monitored  using  one  channel  of 
a  two  channel  recorder  during  the  experiment. 

The  frequency  and  period  of  oscillation  was  determined  by  use  of  an 
acceleroilieter  mounted  on  the  yoke  assembly.  The  signal  was  processed 
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through  a  filter  to  remove  high  frequency,  low  amplitude,  vibrations  and 
was  recorded  on  one  channel  of  a  two-channel  recorder.  Tlie  second  channel 
of  tlie  recorder  was  used  to  record  the  forces  acting  on  the  cylinder. 

The  instantaneous  forces  v/ere  sensed  by  four  piezoelectric  strain  gages 
mounted  on  a  cantilever-beam  type  arrangement,  forming  the  top  of  the 
yoke  assembly.  The  yoke  was  instrumented  so  as  to  measure  both  in-line 
and  transverse  forces,  i  igures  1  tiirouglt  3  show  the  test  section  of  tlie 
tunnel  and  the  cylinder-yoke  assembly. 

System  checkout  and  sensor  calibration  |)receded  the  experimental  runs. 
Several  test  inns  were  made  to  verify  the  proper  operation  of  the  test 
system,  sensors,  processing  and  recording  equipment,  and  to  establish  the 
test  run  procedures  to  be  used.  It  v/as  verified  that  the  mechanical 
system  produced  the  required  periodicity  of  oscillation  for  frequencies 
below  about  four  cycles  per  second  and  amplitudes  of  about  one  inch. 

It  was  verified  that  the  sensors  and  processing  equipment  had  adequate 
sensi tlvity  and  phase  reproduction  to  allow  measurement  of  the  required 
physical  parameters. 

The  first  set  of  experiments  had  as  its  objective,  the  measurement 
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of  tlie  mean  fluid-induced  force  on  the  cylinder,  in  the  direction  of 
the  stream  flow,  (in-line  force)  for  various  frequencies  and  amplitudes 
of  cylinder  oscillation  transverse  to  the  flow.  It  was  verified  by 
experiment  that  there  v/as  no  in-line  component  of  inertial  force. 

The  physical  recording  of  tlie  in-line  force  was  accompl isiied  with 
a  iioneywell  recorder,  i’itot  differential  pressure  was  recorded  slimil- 
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taneously  to  verify  tiie  inaynitutie  and  constancy  of  fluid  velocity  during 
tlie  runs.  Accelerometer  output  was  recorded  on  tlie  two  channel  recorder 
to  provifie  data  on  oscillation  period  and  pliase  angle. 

Subsequently,  tlie  mean  in-line  force,  tlie  flow  velocity,  the  oscil¬ 
lation  period  and  amplitude,  and  tlie  cylinder  diameter  v/ere  recorded 
from  tlie  cliarts.  A  mean  drag  coefficient  C^j.  was  calculated  from  tliis 
data  for  correlation  witli  the  frequency  parameter  D/VT  and  tlie  normalized 
amplitude  A/U  for  the  range  of  Reynolds  numbers  under  consideration. 

The  second  set  of  experiments  had  as  an  objective  the  measurement  of 
the  time-dependent  forces  in  the  direction  of  cylinder  motion,  i.e., 
transverse  to  the  direction  of  flov/.  Since,  in  this  case,  there  was  an 
inertial  force  in  the  direction  of  interest,  it  was  necessary  to  adjust 
the  experimental  procedure  in  order  to  identify  the  instantaneous  value 
of  the  inertial  force.  Separation  of  the  inertial  force  from  the  total 
force  was  required  to  evaluate  the  resistive  force  produced  by  the  fluid 
alone. 

The  total  force  v/as  recorded  simultaneously  with  the  acceleration. 
The  water  level  in  the  tunnel  was  alternately  raised  and  lowered  for  each 
run  in  order  to  allov/  the  recording  of  the  total  "v/et'I  force  and  the 
inertial  "dry"  force  at  a  conwion  amplitude  and  frequency  of  oscillation- 
The  acceleration  trace  provided  the  reference  for  the  time  correlation  of 
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the  wet  and  dry  force  data.  Because  of  the  necessity  to  properly 
correlate  wet  and  dry  phase  information,  extreme  care  was  taken  in 
verifying  the  zeros  for  tlie  acceleration  and  force  recordings. 

From  force  traces  similar  to  Fig.  'I,  the  transverse  force  was 
read  and  recorded  from  hotli  the  v/et  and  dry  run  curves.  Starting  from 
the  zero  acceleration  time  (illustrated  on  Fig.  'I)  the  wet  and  dry 
forces  were  recorded  each  0.01  second,  for  at  least  three  complete  cycles. 

This  data  was  then  punched  on  lUH  cards  for  evaluation. 

I 

The,  second  lest  apparatus  consisted  of  a  recirculating  water  channel. 

The  test  section,  12  inches  v/ide,  18  inches  In'yh  and  36  inches  long,  was 
Closed  on  top  with  a  removahle  plexiglass  plate  to  eliminate  the  free 
surface  effects.  Furthermore,  tv/o  thin  plexiglass  plates  (10  inches  long, 

0  inches  higfi,  and  1/16  inches  thick)  were  mounted  on  the  v/alls  of  the 
test  section,  each  at  a  distance  of  0.0  inches  away  from  the  wall.  The 
cylinder  axis,  at  its  mean  position,  was  at  the  center  of  the  plate.  Tlie 
test  cylinder  v/as  placed  between  these  two  plates.  Tlie  gap  between  the 
plate  and  the  cylinder  end  was  1/32  inches.  The  purpose  of  the  two  plates 
was  to  eliminate  the  effect  of  the  viall  boundary  layers. 

The  drive  system  used  in  connection  with  tlie  recirculating  tunnel 
was  transferred  to  the  water  channel  and  the  test  cylinder  was  connected 
to  it  by  two  thin  leaf  arms.  The  velocity  of  the  fluid  at  the  test  section 
was  regulated  by  adjusting  the  speed  of  the  variable  speed  pump. 

Extensive  velocity  surveys  with  a  hot-film  anemometer  v/ere  made  to 
determine  the  velocity  and  turbulence  distributions.  The  velocity  was 
found  to  be  uniform  in  both  the  vertical  and  horizontal  directions  within 
2%  of  the  mean  velocity.  The  turbulence  level  was  about  0.3%  for  a  mean 
velocity  of  1.5  feet  per  second. 
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As  before,  a  beam-type  accelerometer  was  mounted  on  the  riyid  ann 
connecting  the  twodriving  arms.  The  acceleration  signal  was  electronically 
subtracted  from  the  total  force  signal  by  means  of  a  differencing  unit. 

The  mass  of  the  accelerometer  beam  and  the  sensitivity  of  the  signal 
amplifier  were  adjusted  so  as  to  obtain  zero  net  signal  wljen  the  cylinder 
assembly  was  oscillated  in  air  at  desired  amplitudes  and  frequencies. 

The  natural  frequency  of  the  accelerometer  was  about  80  times  larger  than 
the  largest  vortex  siteddiny  frequency.  Tiie  accelerometer  had  a  logarithmic 
damping  coefficient  of  c  =  0.008.  Tlius,  it  lias  correctly  detected  the 
inertial  force  acting  on  the  system  witlioiit  any  amplification. 

Dynamic  calibration  tests  were  made  by  attacliing  known  masses  to  tlie 
cylinder  and  oscillating  them  at  known  fixed  amplitudes  and  frequencies. 

It  was  found  that  tlie  magnification  factor  of  the  force  transtlucer- 
ampl i fier-recorder  system  vras  nearly  unity,  ami  tlie  natural  frequency 

I 

of  tlie  entire  oscillating  system  was  'IB  llz  in  air  and  '10  llz  in  water. 
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III.  DATA  ANALYSIS 


The  numerical  calculations  as  well  as  measurements  in  time-dependent 
flow  yield  the  resultant  force  as  a  function  of  time  for  a  given  set  of 
numerical  values  of  the  independent  parameters.  Thus,  it  is  not  possible 
v/itliout  a  suitable  hypothesis,  to  express  the  force  both  as  a  function  of 
time  and  remaining  independent  parameters  as  one  ordinarily  v/ould  in  a 
closed  form  solution.  Such  a  working  hypothesis  is  particularly 
necessary  for  the  cable  struiiuiiing  problem  since  Llie  results  are  to  be 
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incorporated  into  tlie  dynamics  of  the  cylinder  motiotj  in  the  form  of  a 
forcing  function.  It  sliould  be  stated  at  the  outset  that  there  is,  at 
present,  no  generally  accepted  hypothesis  to  decompose  the  time  dependent 
force  into  suitable  components. 

Stokes,  in  a  remarkable  paper  on  the  motion  of  pendulums,  showed 
that  the  expression  for  the  force  on  a  sphere  oscillating  in  an  unlimited 
viscous  fluid  consists  of  two  terms,  one  involving  the  acceleration  of 
tlie  sphere  and  the  otlier  the  velocity.  This  analysis  sliows  that  the 
inertia  coefficient  is  modified  because  of  viscosity  and  is  augmented 
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over  the  theoretical  value  valid  for  irrotational  flow.  The  drag 
coefficient  associated  with  the  velocity  is  modified  because  of  accele¬ 
ration,  and  its  value  is  greater  than  it  v/ould  be  if  the  sphere  were 
moving  with  constant  velocity.  In  general,  the  force  experienced  by  a 
bluff  body  at  a  given  time  depends  on  the  entire  history  of  its  accele¬ 
ration  as  well  as  tfie  instantaneous  velocity  and  acceleration.  Thus, 
tl)e  drag  coefficient  in  unsteady  flow  is  not  eciual  to  that  at  the  same 
instantaneous  velocity  in  steady  flow.  Neither  is  the  inertia  coeffi¬ 
cient  equal  to  that  found  for  unseparated  potential  flow.  As  yet,  a 

theoretical  solution  of  the  problem  for  separated  flow  is  difficult  and 
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much  of  the  desired  iiifutiiialioii  must  he  ohtniued  both  experimentally  and 

numerically,  for  example,  Ihrouyli  tlie  use  of  the  discrete  vortex  model. 

In  this  respect,  tlie  experimental  studies  of  Morison  and  his  co-workers 

(1960)  on  the  forces  on  piles  due  to  the  action  of  progressive  waves 

have  shed  considerable  light  on  the  problem.  The  forces  are  divided 

into  two  parts,  one  due  to  the  drag,  as  in  tlie  case  of  flov/  of  constant 

velocity,  and  the  other  due  to  tlie  acceleration  of  the  f’  id.  Tlie 

concept  necessitates  tlie  introduction  of  a  drag  coefficient  and  an 

inertia  coefficient  in  tlie  expression  for  force.  In  particular  if  F 

is  the  force  per  unit  length  experienced  by  tlie  cylinder,  tlien 

« 

F  =  0.5  C^|pl)|lJ|U  I  jipD^/d  .  dl)/dt  (1) 

where  U  and  dU/dt  represent  respectively  tlie  undisturbed  velocity  and  the 
acceleration  of  the  fluid. 

On  the  basis  of  irrotational  flow  around  the  cylinder,  should  be 
equal  to  2  (cylinder  at  rest,  the  fluid  accelerating  ;othenvise  C  =  =  1), 

and  one  may  suppose  tliat  tlie  value  of  should  be  identical  witli  tliat 
applicable  to  a  constant  velocity.  However,  numerous  experiments  siiow 
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tliat  tin’s  is  not  the  case  and  that  C^j  and  show  considerable  variations 
from  those  just  cited  above.  Fven  thougli  no  one  has  suggested  a  better 
alternative,  the , use  of  the  Morison's  equation  gave  rise  to  a  great  deal 
of  discussion  on  what  values  of  tlie  tvio  coefficients  should  be  used. 
Furthermore,  the  importance  of  the  viscosity  effect  has  remained  in  doubt 
since  the  experimental  evidence  published  over  the  said  period  lias  been 
quite  inconclusive. 

The  most  systematic  evaluation  of  the  Fourier-averaged  drag  and 
inertia  coefficients  has  been  made  by  Keulegan  and  Carpenter  (1958) 
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tlirouyli  measurements  on  subineryed  horizontal  cylinders  and  plates  in  the 
node  of  a  standlny  wave,  applyiny  theoretically  derived  values  of  velo¬ 
cities  and  acceleratiojts .  Lxtensive  measurements  by  Sarpkaya  (1976) 
liave  sliown  that  the  dray,  lift,  and  inertia  coefficients  depend  on  botii 
tlie  Reynolds  number  and  tlie  Keulegan-Carpenter  numl)er  (U||T/D)  and  that 
Mori  son's  equation  predicts  remarkably  well  tlie  measured  force  provided 
tliat  the  kinematics  of  tlie  flow  field  is  known  accurately. 

It  is  on  tlie  basis  of  the  foreyoing  that  the  transverse  force 
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acting  on  a  circular  cylinder  undergoing  liarmonic  oscillations  in  the 
transverse  direction  to  a  uniform  flow  is  expressed  in  terms  of  the 
Mori  son's  equation  as 
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where  and  are  tlie  inertia  and  drag  coefficients,  respectively. 
Evidently,  the  normalized  force  as  v/ell  as  the  drag  and  inertia  coeffi¬ 
cients  depend  on  the  parameters  b,„T/n  =  2nA/l),  and  D/VT.  The  first 
parameter  Is  directly  proportional  to  tiie  relative  amplitude,  and  the 
second  parameter  may  be  regarded  as  the  frequency  parameter  D/VT  or  the 
reduced  velocity  Vt/D.  Within  tlie  range  of  Reynolds  numbqrs  considered 
herein  (Re  =  7,000  to  20,000),  in  accordance  witli  the 

specifications  of  the  sponsor,  the  drag  and  inertia  coefficients  are 
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considered  to  be  independent  of  the  Reynolds  number  (Sarpkaya  1976). 

For  the  purpose  under  consideration  it  is  more  advantageous  to 
expand  the  | cos  utjcos  mt  term  in  series  and  to  retain  only  the  first 
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term.  This  procedure  yields. 
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in  which  and  are  given  by  their  Foiiricr-averages  as, 
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In  the  foregoing,  the  inertia  and  drag  coefficients  liave  lieen  denoted 

as  and  C^|  in  order  to  distinguish  them  from  those  corresponding  to 

in-line  oscillations.  The  subscript  'T*  carries  the  meaning  of  "lift” 

or  force  in  the  direction  transverse  to  the  stream,  furthermore, 

and  were  obtained  by  normalizing  F  v/ith  0.5|)l)Uj,j  in  Eqs.  (4)  and  (5). 

—2 

If  the  lift  force  is  noniialized  by  O.bpUV  in  determining  the  dray  and 
inertia  coefficients,  then  Eg.  (2)  reduces  to 

Cl  ^  C,„|,sin  2jrt/T  -  Cj|,cos  2)it/T  (6) 

where  and  are  related  to  and  by 

“  Cji  .[32«A^/D^]/[3V^T^/D^1  ‘  ;  (7) 

and  _  -?  9  ? 

[2ii3A/n]/[vV/0‘']  (8) 

It  should  be  not^d  in  passing  tliat  and  are  identical  to  those 
used  by  llartlen  and  Currie  (1970). 

The  force  acting  on  the  cylinder  in  the  in-line  direction  due  to 
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the  oscillations  in  the  transverse  (iirection  is  expressed  in  terms  of 
a  mean  drag  coefficient,  denoted  by  C^.,  given  as 

^di  **  [Fo>"ce  in  the  in-line  direction]/(0. 5pUV^)  (9) 

The  experimental  data  are  analyzed  according  to  the  Eqs.  ('!)- 
(9)  tlirough  tlie  use  of  afipropriate  computer  programs. 

The  drag  and  inertia  coefficients  obtained  as  described  above  were 
obviously  a  consequence  of  the  use  of  Eq.  (3)  and  tl»e  method  of  Eourier 

I 

averaging.  It  was,  therefore,  necessary  to  compare  tlie  predictions  of 
Eq.  (3)  through  the  use  of  tlie  calculated  coefficients  witli  tliose  measured 
directly.  Evidently,  only  such  a  comparison  could  show  v/hether  the 
decomposition  of  the  instantaneous  force  into  an  in-phase  and  out-of-phase 
component  is  justifiable.  Tlie  calculations  performed  along  these  lines 
have  resulted  in  computer  plots  similar  to  those  shown  in  Figs.  5  through 
9.  Tliere  are  three  curves  on  each  plot:  one  marked  by  COSA  is  the  normalized 
instantaneous  velocity;  one  marked  by  CLTM  is  the  normalized  instantaneous 
transverse  measured  force;  and  finally,  the  one  marked  v/ith  CLTC  is  tlie 
normalized  instantaneous  force  calculated  through  the  use  of  Eq.  (3).  It 
is  evident  that  the  measured  force  is  fairly  well  represented  by  Eq.  (3) 
particularly  for  VT/I)  values  near  critical  synchronization.  Further  away 
from  synchronization  additional  vortex-shedding  frequencies  appear  and 
require  the  consideration  of  higher  harmonics  in  Eq.  (3);.  Such  an  analysis 
has  been  carried  out  but  will  not  be  reported  herein. 


IV.  DRAG  AHO  INERTIA  COEEFICIENTS 


The  results  will  he  discussed  in  t\io  parts.  The  first  will  he  the 
average  in-line  force  acting  on  llte  cylinder  undergoing  forced  periodic 
oscillations  in  the  transverse  direction.  The  second  will  be  the  time- 
dependent  lift  force  and  its  in-phase  and  out-of-pl»ase  cuni|ionents. 

Evidently,  tlie  average  in-line  force  is  coupled  with  secondary 
oscillatjons  due  to  vortex  sliedding.  However,  sucli  oscillations  are 
rather  small  In  both  steady  and  periodic  flows  and  certainly  not  larger 
tlian  about  7%  of  tlie  average  force.  It  Is  for  this  reason  tiiat  ojily  tlie 
average  of  the  in-line  force  acting  on  tlie  oscillating  cylinder  is 
presented  herein. 

Figures  10-12  show  tlie  variation  of  the  normalized  in-line  force 
as  a  function  of  D/VT  for  representative  values  of  A/D.  Each  figure 
represents  the  data  obtained  with  two  velocities,  namely,  V  »  0.04  and 
V  =  1.3  ft/sec.  In  normalized  form,  these  velocities  correspond  to  the 
Reynolds  numbers  Re  =  VD/v  =  7,000  and  Re  =  10,833.  Figure  13  shows 
a  carpet  plot  of  the  mean  drag  coefficient  for  four  representative  values 
of  A/D. 

Evidently,  the  in-line  force  increases  with  A/D  since  the  cylinder, 

undergoing  transverse  oscillations,  presents  a  larger  apparent-projected 

area  to  the  mean  flow.  This,  however,  is  only  part  of  thfe  explanation. 

In  addition,  the  vortex  growth  and  motion  are  affected  by  the  oscillation 
« 

of  the  cylinder  which  in  turn  affect  the  in-line  and  transverse  forces 
acting  on  the  cylinder.  This  is  evidenced  by  the  fact  that  the  in-line 
force  for  a  given  A/0  increases  at  first,  reaches  a  maximum,  and  then 
decreases  as  0/VT  increases.  A  simple  minded  calculation  based  on  the 
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steady  flow  drag  coefficient  for  a  stationary  cylinder  and  tlie  apparent 
projected  area  for  tlie  in-line  force,  which  may  be  written  as 

'C..  =  C  ,  (1  i  2A/r))  (10) 

di  ds 

yields  values  which  are  almost  equal  to  the  inaxiinum  values  given  in 
Fig.  13.  It  sliould  be  noted,  liowever,  that  the  plienomenon  is  far  more 
complex  and  tfiat  such  a  simple  minded  procedure  sliould  not  generally  be 
used,  even  though  the  results  are  surprisingly  good. 

Figure  13  shows  that  tlie  in-line  force  coefficient  reaches  its 

I 

maximum  at  D/VT  between  0.18  and  0.20.  Ordinarily,  tlie  Strouhal  lumitier 
for  a  stationary  cylinder  would  be  0.21  for  the  Reynolds  numbers  cited 
previously,  and  one  would  expect  that  the  forces  acting  on  the  cylinder 
will  undergo  dramatic  changes  as  the  vortex  shedding  frequency  given  by 
the  Strouhal  number  coincides  with  the  frequency  of  the  cylinder 
oscillations.  The  present  results  show  tliat  such  a  synchronization 
takes  place  at  a  frequency  slightly  lower  than  tlie  Strouhal  frequency. 

The  occurrence  of  synchronization  as  well  as  the  increase  of  the  amplitude 
of  oscillations  in  the  force  trace  are  shown  most  dramatically  in  Fig.  14. 

This  figure  was  obtained  by  setting  the  free  stream  velocity  at  0.84  ft/sec 
and  the  A/D  ratio  equal  to  0.5.  Then,  beginning  with  the  case  of  non¬ 
oscillating  cylinder,  the  frequency  of  the  oscillations  was  gradually 
increased  up  to  about  4  llz  and  the  resulting  In-line  force  was  continuously 
recorded.  The  figure  shows  that  the  in-line  force  increases  rapidly  but 
with  very,  little  oscillations  superimposed  on  it.  As  soon  as  the  frequency 
of  oscillations  nears  the  Strouhal  frequency,  the  amplitude  as  well  as 
the  frequency  of  the  force  oscillations  increases.  This  fact  is  worth 

remembering  in  connection  with  the  difficulty  of  obtaining  the  drag  and 
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inertia  coefficients  at  oscillation  frequencies  in  the- vicinity  of  the 
Stronhal  frequency. 

from  an  engineering  viev/point,  tiie  significance  of  the  magnitude 
of  the  in-line  force  is  that  a  cylinder  or  cable  excited  by  the  flow  to 
oscillate  in  the  transverse  direction  may  be  subjected  to  in-line  forces 
several  times  larger  than  those  assume<l  in  its  design,  furthermore,  the 
deflections  caused  by  the  in-line  force  of  sufficiently  flexible  cylinders 
tend  to  couple  with  transverse  oscillations  and  not  only  affect  the 

I 

magnitudp  of  the  transverse  oscillations  but  also  the  path  of  the  cylinder 
motion. 

The  time-dependent  transverse  force  is  described,  as  noted  earlier, 
in  terms  of  a  drag  coefficient  or  C^jj^  and  ar>  inertia  coefficient 
or  [see  Eqs.  (3)  and  (6)J  as  a  function  of  V^=  VT/D.  Figures  15 
through  2^1  show  and  for  representative  values  of  A/l),  (A/l)  =  0.13, 
0.25,  0.50  ,0.75,  and  1.03). 

It  is  seen  from  these  figures  that  important  variations  occur  in 
and  particularly  in  the  vicinity  of  the  Strouhal  frequency  wliere  the 
natural  eddy-shedding  is  both  enhanced  and  correlated  by  the  oscillations. 
The  inertia  coefficient  or  the  normalized  in-phase  component  of  the 
transverse  force  undergoes  a  rapid  drop  as  the  frequency  of  oscillations 
approaches  the  Strouhal  frequency.  In  other  words,  synchronization  or 
lock-in  is  manifested  by  a  rapid  decrease  in  inertial  force  and  a  rapid 
increase  in  the  absolute  value  of  the  drag  force.  Thus,  the  lock-in  is 
a  phase  transformer.  Equally  important  is  the  fact  that  ^ml  readies  a 
value  of  about  2  at  a  value  slightly  under  that  corfesponding  to  the 
perfect  synchronization.  It  has  been  noted  earlier  that  is  equal  to 


unity  for  a  cylinder  oscillating  in  a  quiescent  fluid  and  is  equal  to 

t 

2.0  for  a  flow  oscillating  al)oiit  a  cyliiuler  at  rest.  Thus,  tlie  data 
presented  lierein  show  that  the  net  effect  of  the  cylinder-flow  interaction 
near  sychronization  is  not  unlike  that  of  a  periodic  flow  about  a 
cylinder  at  rest.  In  other  words,  tlie  fluid  becomes  the  oscillator. 

Tlie  data  also  show  that  the  use  of  an  inertia  or  added-inass  coefficient 
equal  to  unity,  as  determined  by  oscillating  the  cylinder  in  a  fluid 
otherwise  at  rest,  is  not  correct  in  modeling  the  vortex- induced 
oscillations.  None  of  the  existing  flow-induced  vibration  models 
predict  the  inertia  coefficietrt  accurately  as  far*  as  tlie  present  data 
are  concerned. 

The  data  for  all  values  of  A/I)  sliow  that  the  peak  value  of 

decreases  from  about  2  to  unity  as  A/D  increases  from  0.75  to  1.03. 

In  other  words,  the  wake  gradually  ceases  i.o  behave  like  an  oscillator 

as  A/f)  approaciies  unity,  for  very  small  values  of  V^,  well  below  the 

synchronization  range,  drops  to  unity,  as  would  be  expected  on 

physical  grounds,  for  values  of  above  that  corresponding  to  the 

perfect  synchronization,  first  drops  sharply  and  then  becomes 

negative.  A  similar  phenomenon  is  observed  in  a  harmonic  flow  about  a 

cylinder  at  rest  (Sarpkaya,  1976),  (or  when  a  cylinder  is  oscillated  in 

a  quiescent  fluid)  for  U^^T/D  values  from  about  0  to  'lO,  depending  on  the 
'■  2 

particular  value  of  D  /uT.  Tf)is  is  a  consequence  of  fractional  vortex 
shedding  as  discussed  by  Sarpkaya  (1976)  and  simply  means  that  tlie  total 
drift  mass  (luring  the  period  of  flow  deceleration  is  larger  than  that 

during  the  period  of  acceleration. 

•*  •  *1 . 

The  dramatic  variation  of  in  tlie  vicinity  of  tlie  synchronization 


region  shows  that  the  fliiid  motion  just  before  the  perfect  synchronization 
is  significantly  different  from  that  just  after  tlie  synchronization.  For 
example,  Griffin  and  Volaw  (1972)  have  observed  that  "vibration  frequencies 
above  and  below  the  natural  frecpiency  f^^  respectively  decrease  and  increase 
tlie  longitudinal  spacing  from  its  value  at  the  natural  frequericy"  and  that 
"both  the  vortex  formation  region  and  the  Itutgi tudinal  spacing  are  influenced 
in  llie  same  way  by  the  frequency,  resulting  itt  atj  expansion  and  contraction 
of  the  wake  at  the  lower  and  higher  frecpiencies." 
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The  drag  coefficient  C^u  or  the  normal  izeil  out-of-phase  compottent 
of  tlie  total  instantaneous  transverse  force  giver>  in  Figs.  15  through  19 
becomes  negative  for  values  in  tht  neighborhood  of  5.  Outside  tin's 
range  tlie  drag  is  mostly  positive,  thus  in  the  opposite  direction  to  the 
motion  of  tlie  cylinder.  Within  the  range  of  values  cited  above,  the 
drag  force  is  in  phase  with  the  direction  of  motion  of  the  cylinder  and 
helps  to  magnify  the  oscillations  rather  than  damp  them  out.  For  this 
reason,  the  region  in  which  is  negative  is  sometimes  referred  to  as 
the  negative  damping  region.  The  fact  of  the  matter  is  that  this  is  not 
damping  in  the  proper  use  of  the  word  but  rather  an  energy  transfer  from 
the  Fluid  to  the  cylinder  via  the  mechanism  of  synchronization.  The 
values  of  at  which  Cjj  or  changes  its  sign  depend  on  A/D.  The 
maximum  absolute  value  of  in  the  synchronization  range  decreases 
rapidly  as  A/D  increases.  Field  studies  have  shown  that  :synchronization 
does  not  occur  for  relative  amplitudes  larger  than  about  unity.  The  trend 
of  the  pi'esent  data  is  in  conformity  with  such  observations. 

Finally,  an  unexpected  and  previously  unknown  observation  in  connection 
with  the  variation  of  will  be  described.  For  normalized  velocities 
in  the  neighborhood  of  4,  the  data  yield  once  again  negative  drag 
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coefficients  fQr  A/D  sinnller  than  about  0.5.  The  occurrence  of  this 
second  region  of  synchroni7ntion  at  higher  fre(|uencies  shows  tliat  there 
is  not  a  single  region  of  lock-in  and  Uiat  there  are  at  least  Iv/o  and 
possibly  more  regions  of  frequency  in  which  synchronization  can  occur. 

The  narrowness  of  the  regions  in  the  second  region  of  synchronization 
makes  it  rather  difficult  to  observe  ti«c  plienomenon.  In  fact,  one  may 
easily  miss  such  a  region  by  simply  not  taking  smaller  increments  in 
frequency.  It  suffices  to  note  that  a  cylinder  may  be  excited  first  at 
frequeniies  near  the  Slrouhal  freipiency  and  then  at  Lite  multiples  of  llie 
Strouhal  frequency.  Ilov/ever,  tlie  largest  energy  transfer  from  tlie  fluid 
to  tlie  cylinder  occurs  in  llie  first  syncronization  region  near  Llie 
Stroulial  frequency. 

The  variations  of  and  with  V^.  for  representative  values  of 
A/D  are  shown  in  Figs.  25  liirough  S^I.  Evidently,  tins  particular  set  of 
data,  obtained  from  and  ,  carry  tlie  same  Information  as  and 

A  careful  examination  of  the  data  shov/s  lliat  in  tiie  synclironization 
region  increases  at  first  gradually  and  tlieii  more  rapidly  with  A/D, 
reaciiing  a  value  of  about  -1.0.  Subsequently,  r.^u^  decreases  rapidly  witli 
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increasing  A/D  and  becomes  positive  near  A/D  1.0.  Similarly,  ^mli 
Increases  with  increasing  A/D  and  reaches  a  mnxiimim  value  of  about  1.65 
at  A/D  =  0.0.  For  larger  values  of  A/D,  decreases  rapidly  as  A/D 
approaches  unity. 

The  values  of  and  corresponding  to  tlie  value  at  the  point 
of  perfect  synchronization  v/ere  entered  into  tlie  computer  and  tlie  equation 
of  motion  was  solved,  as  described  later,  through  the  use  of  a  fourth- 
order  Runga-Kutta  numerical  procedure.  This  melliod  has  enabled  the 

prediction  of  the  dynamic  response  of  cylinders  at  perfect  synchronization. 

,19 
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V.  ANALYSIS 

The  equation  of  motion  for  an  elastically  mounted,  linearly  damped, 
and  periodically  forced  cylinder  may  he  v/ritton  ns  (see  Fig.  35) 


Fig.  35  Elastically  mounted,  linearly-damped  cylinder 

mx  +  cx  >  kx  =  F  =  O.Spj^OLV^Cj^  (11) 

where  m  represents  the  mass  of  the  cylinder,  c  the  linear  material 
damping,  and  k  the  spring  constant.  The  derivatives  of  x  are  taken  with 
respect  to  real  time  t,  i.e..,  d  x/dt  =  x,  etc. 
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equation  (11)  may  be  reduced  to 

X,,  >  ZCX^  ’  x^,  =  p^H^[C^i^sin  Hr  -  (lG/3ii^)  XjC^jjCOS  Hr]  (12) 

in  winch  p^=  Pj-Zp^  =  2Ti^aS^ 

Equation  (12)  nwy  be  solved  exactly  only  for  constant  values  of 
and  Cjjp  There  is  of  course  no  difficulty  in  solving  this  c(|uation 
either  for  assumed  values  of  the  force  coefficients  or  for  theoretically 
or  experimentally  obtained  coefficients.  The  purpose  of  this  report  is 
to  solve  Eq.  (12)  through  the  use  of  the  experimentally  obtained  force- 
transfer  coefficients,  particularly  for  perfect  synciironization,  and 
compare  the  results  v/ith  those  obtained  experimentally  by  others  under 
comparable  conditions. 

Equation  (12)  may  also  be  written  in  terms  of  and  as 


x^  +  ^  x^  =  am^  sin  Hr  -  cos  ill)  (13) 

from  which  one  easily  obtains  the  normalized  relationship  between  the 
amplitude  and  phase  as 

^dh^  "  (M) 

Solving  for  fl,  one  has 

=  (1  -  2c^)  -  t(l  -  2c^)^  •  W  -  1.1^/^  (15) 

where  W  is  given  by 

W  =  (C^„,  +  C,^,,,)/(A/D)^  (16) 


Equations  (13)  and  (14)  constitute  a  set  of  coupled  differential 
equations.  For  a  given  value  of  C^jj^  and  may  be  expressed  in 
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terms  of  A/D  throucjli  the  use  of  the  experi menial  data.  Tlien  the  E(|S.  (13) 
and  (14)  nay  be  solved  toyether  for  given  values  of  'a*  and  C-  To*'  this 
purpose  n  Is  solved  from  Eq.  (15)  for  the  previously  calculated  value  of 
A/I)  and  the  corresponding  values  of  and  Tlien,  Eq.  (13)  is 

integrated  through  llie  use  of  a  fourth  order  Ihjiuja-Kutta  method  for  a 
half  cycle  and  a  new  value  of  A/D  is  calculated.  This  value  is  then 
inserted  into  the  Eq.  (15)  toyether  with  the  corresponding  new  values 
of  and  and  the  entire  procedure  is  repeated  until  a  steady  state 
is  reached.  Tlie  computer  program  wliich  carries  out  the  stated  calculations 
is  given  in  Append! x-A. 

The  purpose  of  the  application  of  the  aliove  procedure  to  self-excited 
oscillations  of  cylinders  was  three  fold:  (a)  to  establish  the  validity 
of  tlie  force-transfer  coefficients  within  the  range  of  Reynolds  numbers 
under  consideration;  (b)  to  investigate  as  to  whether  the  steady-state 
response  of  the  cylinder  is  determined  primarily  by  the  ratio  Sg  =  c/^i 
(see  e.g.,  Skop,  Griffin,  &  Ramlierg,  1976)  or  by  both  C  and  'a';  and 
finally,  (c)  to  examine  the  stability  of  the  differential  equation  v/lien 
the  cylinder  reaches  a  steady-state.  Evidently,  if  it  can  be  demonstrated 
that  the  calculations  (iredict  accurately  the  steady-state  amplitudes  of 
self-excited  oscillations  obtained  under  entirely  different  circumstances, 
then  one  can  conclude  tliat  tlie  force-transfer  coefficients  presented  lierein 
should  form  the  basis  for  comparison  with  those  obtained  from  otlier  numerical 
or  heuristic  models. 

Tlie  foregoing  may  lie  accomplished  only  through  the  use  of  carefully 
conducted  experiments  for  which  tlie  material  damping  ratio  (obtained  in 
vacuo),  the  mass  ratio  'a',  the  critical  frequency  w^,  and  the  steady- 


state  amplitude  are  recorded.  To  the  autlior's  knowledye  only  Griffin 
and  Koopmann  (1977)  were  able  to  obtain  sucb  data  whicli  are  tabulated 
below: 

Damping  ratio  Mass  ratio  Response  parameter  Amplitude 
C  a  S  =  c/a  A/D 

ll 


I 

6.78 

X 

10“'* 

0.62 

xlO"^ 

0.079 

o.n7 

II 

3.31 

X 

10-' 

2.^13 

xio'^ 

O.IO 

0.275 

III 

'  1.62 

X 

io“^ 

5.17 

xKf^ 

0.31 

o.n 

This  data  set  formed  the  basis  of  comparative  calculations 

reported  herein.  Before  proceeding  vntlj  the  calculations,  however,  it 

was  necessary  to  investigate  the  beliaviour  of  the  differential  equation 

and  the  accuracy  of  tlie  numerical  procedure  used,  for  this  purpose, 

C,  and  C  -  v/ere  assumed  to  remain  constant  and  the  differential  equation, 
dl  ml 

[Eqs.  (12)  or  (13)],  v/as  solved  in  closed  form  to  yield 

y^.  =  [sin(-a-<j)/sinp]e  ^  sin(V^-c^  t  »  n)  -  sin(nT-a-<^)  (17) 

in  which  y^  =  x/A  and 

a  =  arc  tan(16ACj^/3n^nC|j^ )  ,  f  =  arc  lan(2cn)/(l-n^) 

and 

n  =  arc tan  [(^?)/[c  t  1)  cot(-«-f)]]  ■  (in) 

Then  the  accuracy  of  tlie  numerical  integration  through  the  use  of  the 
fourth  order  Runga-Kutta  method  has  been  checked  by  comparing  results 
with  those  obtained  from  the  exact  solution  for  various  sets  of  identical 
input  parameters.  In  all  cases  the  solutions  wore  perfectly  identical  in 


both  the  transient  and  steady  states.  Sample  plots  of  the  evolution 
of  the  oscillations  are  shown  in  Fiys.  36  Ihrouyli  30.  Sucli  calculations 
liave  also  sliown  that  a  lowly  damped  system  may  result  In  an  initial 
overshoot  of  the  amplitude  before  it  reaches  a  steady  state  provided 
that  the  force- transfer  coefficients  remain  constant.  A  highly  damped 
syste<ii,  on  tlie  other  hand,  reaches  its  steady-state  gradually  from 
Initial  rest  position.  It  should  be  kept  in  mind,  however,  that  llie 
purpose  of  these  parameteric  calculations  v/as  not  to  investigate  the 

I 

role  played  by  damping  but  ratlier  to  compare  the  predictions  of  the 
exact  solution  witii  those  ol)tained  with  the  numerical  procedure.  In 
a  system  In  which  the  force-transfer  coefficients  vary  witli  amplitude, 
sucli  an  overshoot  may  not  take  place. 

Subsequent  to  the  foregoing  parametric  studies,  computations  were 
carried  out  using  tl>e  data  obtained  by  Griffin  and  Koopmann  (1977)  and 
tile  experimentally  determined  values  of  C^jj^  and  tlirough  tl»e  use  of 
the  numerical  procedure  described  earlier.  The  envelopes  of  the 
translent'and  steady-state  oscillations  for  the  tliree  cases  under 
consideration  are  shown  in  Figs.  39  through  'll.  The  predicted  and 
measured  relative  amplitudes  are  tabulated  below:  ' 


Predicted  A/D 

I 

Measured  A/D,  [Griffin  &  Koopmannl 

S-I 

0.43 

0.47 

0.079 

S-II 

0.235 

0.275 

0.14 

S-III 

0.11 

0.14 

0.31 

Evidently,  the  measured  and  calculated  relative  amplitudes  compare  quite 

i 

well.  The  relatively  small  differences  may  be  attributed  to  the  experimental 
errors  in  the  determination  of  both  A/D  and  the  force-transfer  coefficients. 


/ 


to  the  difficulties  encountered  in  tlie  deteniiinntion  of  the  critical 
damping  factor  in  vacuo,  and  to  tim  variation  in  the  mass  ratio  'a'. 
Furthermore,  it  should  he  noted  that  tlie  relative  amplitudes  reported 
by  Griffin  and  Koopmann  (1977)  are  tlie  maximums  and  tliat  the  mean  values 
of  A/I)  at  synclironization  are  somewhat  lov/er  and  certainly  closer  to 
those  predicted  by  tlie  numerical  analysis. 

It  appears,  at  least  on  tlie  basis  of  the  foregoing,  tliat  tlie  force- 
transfer  coefficients  presented  herein  could  be  used  witli  confidence 

t 

in  the  rjinge  of  Reynolds  mimlters  from  about  fi.riOd  to  25,000  to  predict 
the  dynamic  response  of  elastically  mounted  cylinders  provided  tliat  tlie 
material  damping  and  tlie  mass  ratio  are  given.  Evidently,  additional 
data  are  needed  to  lest  tlie  power  of  prediction  of  the  computer  code 
developed  in  this  investigation  for  all  practically  important  values  of 
Sg  from  about  0.01  to  5.0. 

Next  to  consider  was  the  investigation  of  the  role  played  by  S^, 
known  as  the  stability  or  response  parameter.  Among  others,  Griffin, 

Skop,  and  Ramberg  (1975)  liave  attempted  to  sliow  that  tlie  value  of  A/D 
is  uniquely  determined  liy  the  value  of  a  response  or  stability  parameter 
Sg  v/hich  is  itself  related  to  certain  physical  properties  of  tlie  struiuning 
structure,  i.e.,  =  x,/^.  Vickery  and  Watkins  (1962)  were  the  first  to 

arrive  at  such  a  conclusion  on  the  basis  of  dimensional  analysis  and 
energy-balance  considerations.  Griffin,  Skop,  and  Ramberg  (1975)  obtained 
a  least  squares  fit  to  the  existing  data  (see  Fig.  ^2)  and  arrived  at  the 
following  empirical  formula: 

2A/D  =  1. 29/(1  I  O.'ia  Sp)^*^^ 


(19) 


It  should  be  noted  that  the  2A/D  values  in  I ig.  42  show  large  variations 
(as  much  as  100%)  for  a  given  value  of  Sg  particularly  for  Sg  values 
from  about  0.2  to  1.0,  It  should  also’ be  noted  tliat  the  response 
parameter  used  in  Eq.  (19)  and  in  Fig.  42  was  obtained  by  using  the 
damping  ratios  which  include  both  structural  and  fluid  effects. 

The  equation  of  motion  [Eq.  (13)]  sliows  that  the  response  of  the 
system  is  independently  governed  by  C  and  'a*  and  that  the  said  equation 
cannot  be  expressed  in  terms  of  a  single  response  parameter.  One  can, 
however!  explore  the  reasons  for  the  apparent  correlation  of  tl>e 
experimental  data  with  Sg  (determined  as  noted  above)  through  the  use 
of  Eq.  (14).  For  this  purpose,  let  n  =  1  -  r,  (since  0  is  often  very 
close  to  unity).  Inserting  this  value  in  Eq.  (14),  ignoring  terms  such 
as  e  ,  t;£  *  i.e.,  all  third  and  higher  order  terms,  and  simplifying,  one 
has 

V/n  =  (20) 

in  which  =.  [4,  a 

Evidently,  if  it  is  assumed  that  and  c/a  remain  nearly  constant  then 
and  only  then  one  can  conclude  that  there  is  a  unique  relationship  between 
2A/D  and  Sg.  In  fact,  the  insertion  of  fairly  reasonable  values  of  Cj^, 
and  e/a  into  the  Eq.  (20)  suggest  that  2A/D  may  be  represented  by 

lA/D  =  0.35/[0.12  +  Sg]^^^  (21) 

However,  there  is  no  physical  reason  that  the  two  constants  in  Eq.  (21) 

* 

should  remain  constant  for  all  values  of  the  parameters  cited  above. 

In  fact  the  large  variations  in  2A/D  for  a  given  value  of  Sg  in  Fig.  42 
(particularly  for  Sg  from  about  0.2  to  2.0)  suggest  that  the  lack  of 
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correlation  rather  than  the  experimental  errors  Is  responsible  for  the 
scatter.  Thus,  it  was  deemed  necessary  to  investigate  the  effect  of  Sg 
on  A/D  by  carrying  out  a  parametric  study,  lor  this  purpose,  in  one  set 
of  calculations  C  was  kept  constant  and  'a'  was  varied.  In  another  set 
of  calculations,  'a*  was  kept  constant  and  c  was  varied  so  as  to  arrive 
at  tlie  same  Sg  values.  Ihe  following  table  gives  the  t;  and  'a'  values 
used: 


X 

a 

_ 5/_a _ 

A/I)  (calculated) 

G.DxlO"'’ 

1.36x10“^ 

0.5 

0.23 

M 

6.0  xiO"'* 

1 

0.12 

M 

3.4  xiO"'' 

2 

0.05^ 

'1.3x10"^ 

8.6  xlO"^ 

0.5 

0.36 

8.6x10“^ 

M 

1 

0.18 

1.7x10"^ 

M 

2 

0.11  (unstable) 

As  conjectured  above,  identical  S-  values  did  not  result  in  identical 

U 

A/D  ratios.  However,  surprisingly  enough,  the  A/D  values  for  the  two  se 
of  parameters  did  not  significantly  differ.  Even  though,  additional 
calculations  are  definitely  necessary  to  reach  firmer  conclusions,  one 
can  tentatively  state  that  the  governing  equation  is  not  very  sensitive 
to  the  individual  variations  of  the  damping  and  mass  ratios  provided 
that  Sg  remains  constant.  This  particular  behavior  of  the  differential 
equation  plus  the  use  of  combined  fluid  and  material  damping  may  be 
priiikirily  responsible  for  the  observed  correlation. 

It  is  noted  from  tlie  above  table  that  one  set  of  calculations 
became  unstable  even  though  the  damping  ratio  was  fai»*ly  large  and  an 


error  limit  of  10  ®  was  used  in  the  Runga-KiiLla  integration  sclieme, 

(see  rig.  '13).  It  is  not  clear  wlietlier  the  observed  instability  is 
a  consequence  of  the  numerical  integration,  or  of  the  inherent  instability 
of  the  differential  equation  for  certain  sots  of  the  governing  parameters, 
or  of  the  actual  behavior  of  the  corres|>omling  dynamic  system.  It  will 
be  iiwst  desirable  to  periorm  calculations  with  a  set  of  damping  and 
mass  ratios  for  a  system  which  Iras  exhibited  such  unstable  oscillations 
near  synchronization. 


VI.  CONCLUSIONS 


The  results  presetitecl  l)erein  warrant  the  following  conclusions: 

a.  The  in-line  force  acting  on  a  cylinder  undergoing  haniionic 
oscillations  In  the  transverse  direction  increases  with  A/U.  The  force 
coefficient  reaches  its  tiiaxiiiiuni  at  D/VT  between  0.18  and  0.20; 

b.  The  Fourier-averaged  transverse  force  coefficients  exhibit 
significant  variations  in  the  vicinity  of  the  Strouhal  frequency. 

I 

The  inertia  coefficient  is  larger  than  unity  for  oscillation  frequencies 
larger  than  the  Strouhal  frequency  and  may  be  negative  for  oscillation 
frequencies  lower  than  the  Stroulial  frequency; 

c.  The  variation  of  the  inertia  coefficieiit  with  Vf/D  and  A/l) 
is  not  in  conformity  with  the  predictions  of  the  oscillator  models; 

d.  Tile  use  of  the  experimentally  obtained  drag  and  inertia 

coefficients  for  the  transverse  force  together  with  a  numerical  integrati 
procedure  accurately  predicts  tlie  dynamic  response  of  a  self-excited 
cylinder  in  the  synchronization  region;  ] 

e.  A  parametric  study  of  the  separate  and  combined  effects  of 

t 

material  damping  ratio,  mass  ratio,  and  resfionse  parameter  show  that 
the  maximum  response  of  the  cylinder  is  primarily  governed  by  the 
response  parameter.  Additional  data  and  calculations  are  needed  to 
explore  the  effect  of  tliese  parameters  further; 

f.  It  appears,  on  the  basis  of  somewhat  limited  observations,  that 
the  governing  equation  may  become  unstable  for  certain  combination  of 
the  damping  and  mass  ratios. 
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Fig.  13  Comparison  cun/es:  Mean  in-line  drac  coefficient  versus  0/VT  for 
A/D  =  0.25,  0.50,  0.75,  and  0.3^. 


Fig.  14.  Plienocenological  Deranstraticn:  In-line  force  versus  time  while 
increasing  oscillation  frequency,  (A/Q  =  0.50,  V  =  0.84,  0/VT  * 
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Fig.  ]6  versus  for  A/D  =  0.25 
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Fig.  33  versus  for  A/D  «  0.75  . 
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